Keeping a single qubit alive by experimental dynamic decoupling 
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We demonstrate the use of dynamic decoupling techniques to extend the coherence time of a single 
memory qubit by nearly two orders of magnitude. By extending the Hahn spin-echo technique to 
correct for unknown, arbitrary polynomial variations in the qubit precession frequency, we show 
analytically that the required sequence of 7r-pulses is identical to the Uhrig dynamic decoupling 
(UDD) sequence. We compare UDD and CPMG sequences applied to a single 43 Ca + trapped-ion 
qubit and find that they afford comparable protection in our ambient noise environment. 
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Dynamic decoupling (DD) is a general technique for 
maintaining the phase coherence of a quantum state, with 
particular importance for protecting the quantum infor- 
mation stored in the memory qubits of a quantum com- 
uter The simplest example is the Hahn spin-echo 
, a single 7r-pulse which protects against an arbitrary 
and unknown constant offset in the qubit's precession 
frequency d, 0]. When the state is subject to a time- 
varying offset due to, for example, magnetic field noise, it 
can be protected by a sequence of many 7r-pulses. One of 
these, the Carr-Purcell-Meiboom-Gill (CPMG) sequence, 
is well known in the field of nuclear magnetic resonance 
Q. More recently, other sequences have been investi- 
gated specifically for their dynamic decoupling proper- 
ties, such as Periodic DD, Concatenated DD |(|, random 
decoupling B, composite schemes Q, and local optimi- 
sation (9L lioj: a recent review by Yang, Wang and Liu 



11] 



contains further information and references 

In this paper, we derive a dynamic decoupling sequence 
in a particularly intuitive manner, as an extension to the 
spin-echo 0. We prove that with n pulses, the sequence 
can cancel out all the dephasing that would be caused by 





FIG. 1: Illustration of the effect of dynamic decoupling se- 
quences on the acquired phase, as calculated by ([2]). a): Fre- 
quency offset 8 varies linearly with time, and the phase shift 
can be completely corrected by a sequence with two 7r-pulses 
(for n = 2, CPMG and UDD are identical), b): 5 varies as 
a (arbitrary, unknown) quartic polynomial, and is perfectly 
corrected by n = 5 pulse UDD. 
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the frequency varying as an (n — l)th order polynomial 
function of time, without knowledge of the polynomial 
coefficients. This sequence is identical to the Uhrig Dy- 
namic Decoupling (UDD) sequence 12|, |l3j, which was 
originally derived by considering the interaction of a spin 
qubit with a bosonic bath. We implement the sequence 
on a sins le 43 Ca+ ion, demonstrating that the coherence 
time of this qubit is significantly increased, and compare 
it with the CPMG sequence. 

Suppose an arbitrary qubit state is prepared at time 0, 
and we want to recover it at time r. The pulse sequence 
is a series of (assumed ideal and instantaneous) 7r-pulses 
at times ot\T, a-2j, . . . , a n r, where the on are to be found. 
We have remarked that a single Hahn spin-echo will cor- 
rect for a constant frequency offset. If the offset varies 
linearly with time, we can correct the phase error with 
two 7r-pulses at t = \ and f, where t — time/r (Figure 
QJ,). To generalise further, postulate that n pulses suf- 
fice to correct for a frequency variation 8(t) that is an 
(ri — I)th-ordcr polynomial in time (Figure [TJ)) : 



6(t) = Po + pit + p 2 t 2 + ■■■+ p^xf 



(1) 



The phase error <j) eTT is given by integrating 5(t) over time. 
But each 7r-pulse reverses the direction of the qubit's pre- 
cession, so between pulses i and i + 1, if i is odd, we mul- 
tiply the acquired phase by (—1). The resulting integral 
is thus 

</w =Er =0 (-i)7 Q Q / +i ^)di 



oti+i 



(2) 



where ao = and a n +i — 1. Collecting terms for each 
polynomial coefficient pf 
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We require <j) err to be for any choice of the pj , and so 
we obtain a set of n simultaneous equations for the ctj 



(-1)*- 2j2(-iyc4=0 Vj = l,2, 
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These are solved by 



2 n + 1 



(5) 



which can be proved directly by substituting into (fj| 
and applying a series of trigonometric identities [14j | . 

The sequence is independent of r; however in prac- 
tice the frequency offset S(t) is only approximated by a 
polynomial, and as r increases we need more polynomial 
terms (and hence more 7r-pulses) for the approximation 
to be valid. 

This sequence was previously and independently dis- 
covered by Uhrig [12|, |l3j, by considering the spectral 
properties of a qubit coupled to a bath of bosons that 
cause decoherence. The echo sequence was treated as a 
filter in frequency space. Uhrig demanded that the first n 
derivatives of the filter function vanish at zero frequency, 
because this gives the strongest suppression of the noise 
at low frequencies, and this condition leads to the simul- 
taneous equations ^ and hence the sequence JS]). Lee, 
Witzel and Das Sarma have shown [l5| that this sequence 
is optimal for any dephasing Hamiltonian, where "opti- 
mal" means that it is the sequence that maximises the 
qubit fidelity in the small r limit, for a given number of 
pulses [31(| . While this paper was in preparation, Hall et 
al. have independently published a derivation equivalent 
to ours 1161. 



In a 1988 paper [17|, Keller and Wehrli suggest using 
a theoretical procedure similar to ours, to cancel the ef- 
fects of successive polynomial orders of fluid flow in MM. 
However, this "gradient moment nulling" allows 6(t) to 
be controllably scaled by the experimenter; Keller and 
Wehrli concentrate on this parameter rather than pulse 
timing and so do not find the UDD sequence. 

The first experimental tests of UDD were by Biercuk et 
al. , who applied a variety of dynamic decoupling schemes 
to ensembles of ~ 1000 9 Be + ions in a Penning trap fljj . 
Dynamic decoupling was demonstrated in a solid by Du 
et al. (using electron paramagnetic resonance of ensem- 
bles of unpaired carbon valence electrons in irradiated 
malonic acid crystals) [l9[ , and in a dense atomic gas by 
Sagi, Almog and Davidson (~ 10 6 87 Rb atoms in a dipole 
trap) [20(. And recently, Ryan, Hodges and Cory imple- 
mented sequences using single nitrogen vacancy centres 
in diamond 2l| . 

We have applied dynamic decoupling to a single 43 Ca + 
trapped- ion qubit, held in a radio- frequency Paul trap 
[22j. The qubit is stored in two hyperfine states in the 
ground level, \\) = 4S^ 4 and |t> = 4S"}+ 3 (where the su- 
perscripts indicate the quantum numbers F,Mp); these 
states are separated by a 3.2 GHz Ml transition. The 
transition's sensitivity to the external magnetic field is 
2.45 MHz G" 1 at low field; we apply a field of 2.2 G to 
define a quantization axis and to increase the ion's flu- 
orescence rate (by destabilising dark states 23J). Rabi 



oscillations are driven on the qubit transition at Rabi fre- 
quency 2tt x 18 kHz, using microwaves. These are gener- 
ated using a versatile synthesizer, amplified with a solid- 
state amplifier (to » 750 mW) and broadcast inside the 
vacuum chamber using a trap electrode as the antenna. 
To improve the fidelity of the dynamic decoupling tt- 
pulses we apply a small 50 Hz signal, synchronized with 
the AC line, to a magnetic field coil which cancels the 
dominant component of the magnetic field fluctuations 
experienced by the ion; the remaining noise has ampli- 
tude up to ±3 kHz. Each experimental sequence is also 
line-triggered. 

Each experiment (Figure [2^) starts with the ion op- 
tically pumped into state A decoupling sequence 
is tested by sandwiching it between two ^-pulses. The 
second pulse has a phase offset relative to the first; 
scanning this phase leads to Ramsey fringes. Any loss 
of phase coherence in the Ramsey gap leads to fringes 
of reduced contrast, so generally the contrast falls as the 
gap is made longer. We aim to show that this fall be- 
comes slower when dynamic decoupling is used. Finally 
the qubit state is measured by electron shelving and flu- 
orescence detection, with accuracy up to 99.8% [zj ]. 

The sequence is repeated 200 times for each value of 
<fi, which is typically scanned from —450° to +450° in 20 
steps resulting in the measured state varying sinusoidally 
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FIG. 2: a) Experimental sequence: the qubit coherence after 
time r is measured by a Ramsey experiment in which the 
phase <f> of the second -|-pulse is scanned relative to that of the 
first. An n = 4 pulse UDD dynamic decoupling sequence is 
shown in the example (the ~ -pulse and 7r-pulse durations are 
exaggerated for clarity), b) Data (with shot noise error bars) 
and fitted Ramsey fringes for n — 3 and n = 6 pulse UDD 
sequences, both at r = 7 ms. The dramatic improvement in 
the qubit coherence given by the 6-pulse sequence is clear. 
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FIG. 3: Comparing UDD sequences with from to 20 ir- 
pulses. The error bars indicate the uncertainty in each fringe 
contrast fit (estimated by a bootstrap algorithm). The curves 
show the theoretically predicted contrast for each sequence, 
using ([7J and the fitted S(oj). Multiple data points with the 
same r have been combined for clarity. 
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FIG. 4: Comparing UDD and CPMG sequences for six and 
twenty 7r-pulses. Solid symbols (error bars omitted for clarity) 
and dotted lines are the results and fits from UDD sequences 
as shown in Figure [3] Hollow symbols represent CPMG se- 
quences, with solid lines the theoretical prediction using the 
same fitted noise spectrum S(lu) as for the UDD sequences. 
Data points with the same r have been combined for clarity. 



with 4>. A sine curve is fitted to the data to measure 
the contrast; example data is shown in figure [5}). Typ- 
ically 10-20 such runs are taken for a given decoupling 
sequence, with t chosen to be different for each, and with 
the decoupling pulse timings being scaled accordingly. 

Figure [3] shows the results for different numbers of 
7r-pulses. With no dynamic decoupling 7r-pulses, the 
fringe contrast drops to 1/e of its initial value in a time 
t c = 0.51(5) ms; with a 20-pulse UDD sequence, this 
time is extended to r c = 33(1) ms. We also compared 
UDD and CPMG (equally spaced 7r-pulses, at times 
oti = (i — 1/2) jn for i = 1 . . . n) sequences, with results 
shown in figure [4] It can be seen that, in our noise envi- 
ronment, UDD performs no better than CPMG; indeed, 
CPMG is slightly better, extending the coherence time 
to t c = 37(1) ms for a 20-pulse sequence, an increase over 
the unprotected qubit by a factor « 73, or 1.9 orders of 
magnitude. The similar performance of UDD and CPMG 
is expected if the noise spectrum extends to high frequen- 
cies; UDD would be superior if the noise spectrum had a 
sharp high-frequency cutoff 18, 25|. We also performed 
experiments both with and without a 90° phase shift on 
the UDD 7r- pulses [32|, which is equivalent to testing the 
dynamic decoupling for two different qubit states on the 
equator of the Bloch sphere; there was no significant dif- 
ference between the results. 

To fit the data in figures [3] and |4l we perform a simu- 
lation based on the filter function formalism of Cywinski 
et al. 26]. Suppose that the noise power spectrum is 
given by S(uj). We multiply the noise spectrum by the 
pulse-sequence's filter function F(ujt) (F also depends on 
the number and finite duration of the 7r-pulses 18]), and 



calculate the integral over angular frequency w: 

S(u)F(ujt) 



x(t) 



-duj 



The qubit coherence C(t) is then given by [26 
C(t) = Ne- x( $ 



(6) 



(7) 



where N is a normalization constant that accounts for 
effects such as imperfections in the 7r-pulses themselves. 
In our experiment, C(t) is the contrast of the Ramsey 
fringes. 

The noise spectrum of the magnetic field measured out- 
side the ion trap vacuum system did not give a good 
fit to the data when used to calculate C(t), presumably 
because it differs too greatly from the noise at the po- 
sition of the ion. However, we can reverse the process; 
the dynamically-decoupled ion acts as a spectrometer to 

H, [Ijj. 



We model the 



measure the field fluctuations 
noise spectrum S(oj) by a piecewise cubic spline in log- 
log space, use it to calculate C(t), and find the spectrum 
which gives the best fit to the experimental data; the fit 
attempts to match all our UDD and CPMG data with the 
same S (oj) (though each data set is allowed its own fitted 
normalization constant N). The calculated contrast C(t) 
is not very sensitive to the detailed shape of S(u>), but 
the procedure does yield a noise spectrum which is close 
to a power law for 100 Hz < (cj/2tt) < 100 kHz, with 
S(ui) oc ui~ 5±1 . This is consistent with the S(w) oc uj^ 4 
spectrum measured by Biercuk et al. inside a supercon- 
ducting solenoid [§] . The curves in figures [3] and |4] show 
the calculated C (t) using this noise spectrum, and fit the 
experimental data reasonably well. 
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FIG. 5: The 1/e coherence time, as measured by ((7}, is plotted 
against the number of 7r-pulses for UDD and CPMG data. A 
straight line fit to all the data is shown. When the Ramsey 
delay r is equal to the coherence time, the finite-length tt- 
pulses occupy a total duration < 0.02r in all cases. 



The fitted 1/e coherence times are shown in figure [5] 
The data is matched well by a straight line, similar to 
the observations of Ryan, Hodges and Cory 21 1. 

It is clear from figure [3] that although the UDD se- 
quence significantly extends the coherence time of the 
qubit, the coherence at short time is actually degraded 
due to imperfections in the 7r-pulses which are more sig- 
nificant the larger the number of pulses used. We esti- 
mate the typical 7r-pulse fidelity (based on the fits ex- 
trapolated to t = 0) to be 98.7%. This fidelity could be 
improved significantly by increasing the Rabi frequency 
so that it is well above the amplitude 6(t) of the dominant 
noise sources, for example by driving the qubit transition 
with near-field microwaves from electrodes much closer to 



the ion, as proposed in [27 1. 

In conclusion, we have shown that extending the Hahn 
spin-echo to correct for frequency offsets which vary poly- 
nomially in time yields the Uhrig dynamic decoupling se- 
quence, and that applying this sequence (or the CPMG 
sequence) to a single physical qubit stored in a trapped 
43 Ca + ion increases the coherence time by nearly two or- 
ders of magnitude, to r c « 35 ms. In order to demon- 
strate the increase in coherence time, we chose qubit 
states in the S 1/2 manifold which had the greatest sensi- 
tivity to magnetic field fluctuations. For a qubit stored 
in the magnetic field-insensitive "clock" states (4S^° and 
4S*y° 2 ) we have previously measured a coherence time 
T 2 = 1.2(2) s since this was also limited by magnetic 
field noise, it should be possible to extend the coherence 
time of such a qubit to several minutes using dynamic 
decoupling techniques, at which point it becomes practi- 
cally difficult to measure using a single qubit. The mem- 
ory qubit coherence time would then exceed the typical 
timescale for trapped-ion quantum logic gates (~ 20 /is 



[28|) by many orders of magnitude, an essential prereq- 
uisite for implementing fault-tolerant quantum computa- 
tion. 
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